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Abstract 

We carry out a quantum analysis of a dc SQUID mechanical displacement detector, comprising a 
SQUID with mechanically compliant loop segment, which is embedded in a microwave transmission 
line resonator. The SQUID is approximated as a nonlinear, current dependent inductance, induc- 
ing an external flux tunable, nonlinear Duffing self-interaction term in the microwave resonator 
mode equation. Motion of the compliant SQUID loop segment is transduced inductively through 
changes in the external flux threading SQUID loop, giving a pondcromotive, radiation pressure 
type coupling between the microwave and mechanical resonator modes. Expressions are derived 
for the detector signal response and noise, and it is found that a soft-spring Duffing self-interaction 
enables a closer approach to the displacement detection standard quantum limit, as well as cooling 
closer to the ground state. 

PACS numbers: 85.25.Dq; 85.85.+j; 03.65.Ta 



I. INTRODUCTION 



Recently there has been interest in exploiting the nonlinear dynamics of nanoelectrome- 
chanical systems (NEMS) for amplification .'^'^ The use of nonlinear mechanical resonators 
to some extent parallels investigations with systems comprising purely electronic degrees 
of freedom, such as nonlinear superconducting devices incorporating Josephson Junctions 
(JJ)™For example, it was shown that the bistable response of an RF-driven JJ can be 
employed as a low noise, high-sensitivity amplifier for superconducting qubitsP A similar 
setup consisting of a JJ embedded in a microwave cavity was used to measure the states 
of a quantronium qubit,'^ where the relevant cavity mode was found to obey the Duffing 
oscillator equationP 

One area of nanomechanics that has yet to fully explore the possibility of exploiting non- 
linearities for sensitive detection involves setups in which a nanomechanical resonator couples 
either capacitivelj^^'^^'^ or induct ivelj^^^'^ to a superconducting microwave transmission line 
resonator, combining elements from both the above-described NEMS and superconducting 
systems. Such setups are in some sense the solid-state analogues of optomechanical sys- 
tems, which ponderomotively couple a movable mirror to the optical field inside a cavity 
using radiation pressureJ ^^ * ^^ * ^^ * ^^ * ^^ * ^^ * ^" * ^^ ^ In both areas, the focus has primarily been on 
operating in the regime where the cavity and resonator behave to a good approximation as 
harmonic oscillators interacting via their mutual ponderomotive coupling. However, in the 
case of microwave cavities, introducing an embedded JJ,'Slor simply driving the cavity close to 
the superconducting critical temperature,'^ results in the breakdown of the harmonic mode 
approximation; nonlinear dynamical behavior of the cavity must be taken into account. Fur- 
thermore, the ponderomotive coupling term between the microwave or optical cavity mode 
and mechanical mode is by itself capable of inducing strong, effective nonlinearities in the 
respective mode equations. In optical systems, such nonlinearities can manifest themselves 
in the appearance of a bistable (or even multistable) region for the movable mirror.'^^'^I! gy 
restricting ourselves to linear microwave cavities, we are overlooking a range of nonlinear 
phenomena that might enable a closer approach to quantum-limited detection, as well as 
cooling of the mechanical oscillator closer to its ground state. As an illustration, consider 
the phase sensitive Josephson parametric amplifier ,'2SI2Z11SI29] -^j-^ich exploits the nonlinear ef- 
fective inductance of the JJ to perform (in principle) noiseless amplification and quantum 
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squeezing of the respective complimentary quadrature amplitudes of the signal oscillator. 

In this paper, we will go beyond the usually considered ponderomotively-coupled two 
oscillator system to include a Duffing nonlinearity in the microwave cavity mode equations. 
The closed system model Hamiltonian describing the nonlinear microwave-coupled mechan- 



ical oscillators is given by Eq. (26). The nonlinear microwave mode is externally driven 
with a pump frequency uOp that can be detuned from the transmission line mode frequency 
ujt- Our investigation will focus on the nonlinear amplifier created by embedding a dc- 
SQUID displacement detector into a superconducting microwave transmission line.!^ This 
has the advantage of significantly larger coupling strengths-^ as compared with existing ge- 
ometrical coupling schemesJ ^ * ^" * ^^ ' The displacement detector comprises a SQUID with one 
segment consisting of a doubly-clamped mechanical resonator as shown in Fig. [T} The net 
fiux, and therefore circulating current, is modulated by the mechanical motion, providing 
displacement transduction. The capacitively-coupled pump/probe feedline both drives and 
provides readout of the relevant transmission line resonator mode amplitude (or phase). We 
will assume transmission line losses are predominantly due to coupling with the feedline, 
and that the pump drive is coherent. The main irreducible noise source is therefore mi- 
crowave photon shot noise from the drive that acts back on the mechanical oscillator via 
the intermediate nonlinear microwave resonator and SQUID. Environmental infiuences on 
the mechanical oscillator other than that due to the SQUID detector are simply modelled 
as a free oscillator thermal bath. By operating the amplifier well below the superconducting 
critical temperature, and with transmission line currents less than the SQUID JJ's critical 
current threshold, resistive tunneling of electrons and the associated noise is a negligible 
contribution. Similar setups involving JJ elements have been considered previously.'^^HSEIEa 
With JJ plasma frequencies assumed to be larger than both the mechanical and trans- 
mission line fundamental mode frequencies, the SQUID can be considered as a passive, 
effective inductance element that depends on both the externally applied fiux and drive 
current. The effective inductance can therefore be freely tuned by varying these external 
parameters. Previously, we considered only the lowest, zeroth order expansion of the in- 
ductance with respect to the current entering (or exiting) the SQUID,'^ yielding the usual 
ponderomotively-coupled double harmonic oscillator system. In this companion paper, we 
include the next leading, quadratic order term, resulting in a nonlinear current dependent 
inductance. Provided that the current magnitude is small as compared with the JJ's critical 
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current, neglecting higher order terms should not introduce significant errors. The nonlinear 
inductance induces an effective Duffing (i.e., cubic) self-interaction term in the microwave 
mode equations of motion. The results presented here apply to a broad class of bosonic 
detector, which includes optomechanical amplifiers with nonlinear cavitie^ that are de- 



scribable by Hamiltonian (26). A related analysis of quantum noise in a Duffing oscillator 
amplifier is given in Ref . |3ll 

The paper is organized as follows. In Sec. |TT] we first derive the truncated Hamilto- 



nian (26) that describes the closed system dynamics of the coupled cavity and mechanical 
resonator fundamental modes. We then derive the quantum Langevin equations of motion 
that describe the open system dynamics in the presence of the pump/probe line and me- 



chanical oscillator's external environment. In Sec. Ill we find expressions for the detector 
signal response and noise using a semiclassical treatment of the detector's linear response 



to the external noise input signal driving the mechanical oscillator. Section IV determines 
the critical drive current for the onset of bistability (not to be confused with the JJ critical 
current). Sections |V] and VI discuss the effects of the microwave mode Duffing nonlinear- 
ity on mechanical mode displacement detection and cooling, respectively, giving illustrative 
examples assuming achievable device parameters. Section VII| briefiy concludes, while the 
more technical aspects of the signal and noise term derivations are relegated to the appendix. 



II. EQUATIONS OF MOTION 
A. Closed System Hamiltonian 

The displacement detector scheme is shown in Fig. [T| The device consists of a stripline 
resonator (transmission line) of length I bisected by a SQUID. The transmission line is char- 
acterized by an inductance and capacitance per unit length Lt and Ct respectively. The 
SQUID comprises two identical Josephson junctions with critical current Ic and capacitance 
Cj. A segment of the SQUID loop is mechanically compliant, forming a doubly clamped 
resonator of length /qsc- We only take into account mechanical fundamental mode displace- 
ments in the plane of the loop and assume that the resonator can be modeled effectively as a 
harmonic oscillator with the y coordinate giving the centre of mass displacement. The mag- 
netic ffux threading the loop is given by $cxt(l/) = $cxt(0) + XBe^tLscV, where $ext(0) = $ext 



4 



is the flux with the mechanical oscillator flxed at y = 0, Bcxt is the externally applied fleld in 
the vicinity of the resonator, and A is a geometrical factor that corrects for the non-uniform 
displacement of the oscillator along its length. The coupling between mechanical oscillator 
and external heat bath is characterized by the oscillator amplitude damping rate 'jbm, while 
the pump-probe line-transmission line coupling is characterized by the transmission line am- 
plitude damping rate ■jpT- In what follows, we will assume weak couplings (i.e., large quality 
factors for the transmission line and mechanical oscillator) and also that the dominant dis- 
sipation mechanism for the transmission line is due to its coupling to the pump-probe line, 

IpT- 




-l/2< i ►VS 



FIG. 1: Layout of the dc SQUID displacement detector. The dimensions of the SQUID have been 
enlarged relative to the transmission line to show the key characteristics employed in the analysis. 

In analyzing the SQUID dynamics, an appropriate choice of variables is 7-1- = (0i ±02)/2, 
where 0i and 02 are the gauge invariant phases across the Josephson junctions,'^ while for 
the transmission line we choose the phase fleld (j){x,t). The transmission line current and 
voltage are described in terms of t) using the standard telegraphic relations: 



Vi., ^ (2) 

where $o = h/{2e) is the flux quantum. Assuming that the SQUID can be lumped at the 
midpoint x = of the transmission line, the equations of motion for the closed system 
comprising the SQUID, transmission line and mechanical oscillator are given bj^ 

d'^S 1 d'^S 



^/l- + cos(7+) sin(7_) + 2/?^^ 



7_ — TT n + 







(3) 

(4) 



-2- ^ ■ / ^ f , , "^0 '90(0, t) 
uj, 7+ + sm(7+) cos(7-) + ^^^^^ = 0, (5) 



and 



my + mujl^y - ^XBe^JoscI- = 0, (6) 



where uj = ^JlT^Id (Cj$o) is the plasma frequency of the Josephson junctions, /3l = 
l-nLl^j^Q is a dimensionless parameter with L the SQUID loop self-inductance and 1^ the 
Josephson junction critical current, and where n takes on integer values arising from the re- 
quirement that the phase around the loop be single- valued. Equation ([3]) is the wave equation 
for the transmission line, equation Q describes the current circling the loop, which depends 
on the external flux and oscillator position, equation (|5]) describes the average current in the 
loop, and Eq. ([6]) is Newton's second law for the mechanical oscillator with Lorentz force 
acting on the oscillator. The current and voltage across the SQUID must also obey the 
boundary conditions 

g0(±//2,t) _^ c}0(O-,t) _ c}0(O+,t) 
. L gV(0,t) _ g0(O-,t) g0(O+,t) 

^+ 4Lt dtdx dt dt ' ^ ^ 

Using Eqs. Q-Q, we shall now derive approximate equations of motion describing a 
single mode of the transmission line interacting with the mechanical oscillator, where the 
form of the interaction between the two oscillators is governed by the SQUID parameters 
and boundary conditions. Assume that the following conditions are satisfied: (a) uj 3> 
ujt > i^m- (b) (3l < 1. (c) |/(0,t)//c| < 1. (d) |A5ext^oscy/$o| < 1- Condition (a) states 
that the SQUID plasma frequency is much larger than the transmission line mode frequency 
of interest, cut, and consequently we shall ignore the SQUID inertia terms in Q and 



Condition (b) allows us to neglect the SQUID loop self inductance and, together with (a), 
eliminate 7-1- from the equations by expressing them in terms of the transmission line and 
oscillator coordinates as series expansions in P^. Conditions (c) and (d) allow us to expand 
the above equations in the transmission line current l{0,t) = I{t) at x = and in the 
oscillator displacement y. Keeping terms to first order in y and to leading, second order in 
/, Eq. (|6| for the mechanical oscillator becomes approximately 

my + muly - LoiP/2 = 0. (9) 



(10) 



The voltage boundary condition (|8j) can be expressed as 

d $0 \ d<p{o-,t) g0(O+,t) 

^[^a>l/K] = ^[^^ 

where L{I,y) is the effective inductance, which expanded to second order in I takes the 
form 

L{I, y) = Loo + [I/Icf + L,,y, (11) 
where the Lij coefficients are defined as 

$0 

^00 = — — sec (7r<l>ext/$o) (12) 
L20 = -^sec3(vr$ext/$o) (13) 
= ^^^1^ sec (7r$e.t/$o) tan (vr^ext/^o) • (14) 
Note that we have neglected the Py term in (11), restricting ourselves to the leading order 



coupling only between the transmission line and mechanical oscillator, as already stated. The 
above equations differ from those of the prequeP^ through the inclusion of the nonlinear, 
leading order current-dependent contribution [L2q{I / If)"^] to the effective inductance L(/, y). 



The nonlinear voltage boundary condition (jlOj) with inductance given by Eq. (11) gener- 
ates frequency tripling harmonics of the transmission line resonator mode. Omitting for the 
time being the mechanical oscillator degree of freedom, a trial perturbative mode solution to 
the wave equation ^ that includes the leading harmonic and solves the current boundary 
conditions ([T]) is the following: 



0(x,t) 



+A cos{ujt + if) cos [k{x - 1/2)] + aA^ cos(3u;i(: + 3ip) cos [3k{x - 1/2)] ; x > 

-A cos{iut + ip) cos [k{x + 1/2)] - aA^ cos(3a;t + 89?) cos [?,k{x + 1/2)] ; a; < 0, 

(15) 



where k = fc^"-* + k^^'' and cu = \k\l \fLrpCT- The coefficients a, k^^^ and k^^^ are determined 



by substituting Eq. (15) into the vohage boundary condition ViSM and solving perturbatively 



to order A? ^ with k^^^ scahng as . We obtain: a = —1/48, 

(fc(°)Z/2) tan (A;(°)Z/2) = C~' (16) 

and 

k^l = -1(3^2(^(0)^)3 gj^2 (^(0);/2) , (17) 

where 

Considering the transmission hne phase field at the location x = —1/2, where the field is 
pumped and probed (see Fig. [T]), the perturbative solution (15) can be obtained from the 
following single mode equation for (f){—l/2,t) = (pit): 

^ + 40 + ^4 (l - 18C^ [(A:(°)//2) sin (fc(°)//2)]^) 0^ 

+ 1 (l - 2C^ [(A:(°)//2) sin (fc(°)//2)]^) ^ = 0, (19) 

where ujt = \k'^^^\/y/Lj'CT- The awkward nonlinear term 0^ be eliminated by redefining 
the phase mode coordinate as = -0(1 + Tip'^), provided |r|0^ -C 1, where 

r = (l - 2C' [{k^'h/2) sin (A;(°)//2)]') . (20) 



The mode equation (19) in terms of the redefined phase coordinate ijj then becomes 



+ 4^ - l^C' [(A;(°)//2) sin (A;(o)V2)]'^' = 0. (21) 

Thus, embedding a SQUID in a microwave transmission line induces a cubic nonlinearity 
in the effective single mode equations (under the conditions of small currents as compared 
with the Josephson junction critical current), resulting in the familiar (undamped) Duffing 
oscillator. 

We now restore the mechanical degree of freedom y{t) by assuming that for small and 
slow displacements [conditions (a) and (d) above], the interaction with ip can be obtained 



by expanding lut [through its dependence on $ext(2/)] to first order in y in Eq. (21 ) to obtain 



i^ + uj^j.^ - -4C'[(fc^°^V2)sin(A;W//2)]'7/>3 
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tan (7r$ext/$o) sec (7r$ext/$o) uj^'ip- 



(22) 



($o/7r) 47rLr//c 

Equation ^ for the mechanical oscillator, together with Eq. (22) for the phase coordinate, 
follow from the Lagrangian: 



J - ^f"'°f , tan (^*„/*„) sec (^*,.,/4.„) 



uj^ip (t) 



-4c' [(A:(°^V2) sin(A;(°)//2)]'^^ 



(23) 



Introduce the phase momentum coordinate = dC/dijj = m^ip and raising (lowering) 
operators 

1 



1 



m^UTip =F «Pv 



(mcuml/ =F iPy) 



(24) 
(25) 



satisfying the usual commutation relations, where the effective phase mass is = 
|C2./($o/27r) sin2(A;o//2). In terms of the raising (lowering) operators, the Hamiltonian 
operator is 

E = hwro^aT + ^huxKaia:^ + at)^ + fi^mati^'m + ]^h^TKTm{a'T + 0'tY{o.'L + «m), (26) 

where, for notational convenience, hats on the operators and the minus superscript on the 
lowering operators will be suppressed from now on. The parameter characterizing the 
strength of the interaction between the transmission line mode and mechanical oscillator 
mode is 

$0 



K' 



Tm 



tan (7r$cxt/$o) sec (7r$ext/$o) 



(27) 



where A^p = ^Jfij (2mujjn) is the zero-point displacement uncertainty. The parameter 
characterizing the strength of the Duffing nonlinear term takes the form 



(2e)V(2CT0 
huJr 



(28) 



which has been written in such a way as to make clear its various dependencies. In particular, 
Kd depends essentially on the cube of the ratio of the linear SQUID effective inductance 
Lqq to transmission line inductance L^Z, as well as on the ratio of the single Cooper pair 



charging energy to the microwave mode photon energy of the transmission hne. Since the 
strength and sign of the hnear SQUID inductance depends on the external flux $cxt [see 



Eq. (12)], it is possible to vary the strength as well as the sign of the Duffing constant by 



tuning the external flux either side of $o/2. Thus, we can have either spring hardening or 
spring softening of the transmission line oscillator mode. Previously this flux tunability was 
observed in the readout of a persistent current qubitPNote, however, that the perturbative 



approximations that go into deriving the above Hamiltonian (26) do not allow too close an 



approach to the singular half-integer flux quantum point. In particular, the validity of the 
expansions in It and 13 l properly require the following conditions to hold: 

/ 



sec (7r<l)ext/$c 



< 1 



|/?Lsec(7r<l>cxt/<l'o)| < 1- 



(29) 
(30) 



As already noted, Eq. (26) without the Duffing nonlinearity coincides with the Hamil- 



tonian commonly used to describe the single mode of an optical cavity interacting with a 
mechanical mirror via the radiation pressure. However, we have just seen that embedding a 
SQUID within a microwave transmission line cavity induces a tunable Duffing self-interaction 
term as well; it is not so easy to achieve a similar, tunable nonlinearity in the optical cavity 
counterpart. 



B. Open System Dynamics 

Up until now we have considered the transmission line, SQUID and mechanical oscillator 
as an isolated system. It is straightforward to couple the transmission line to an external 
pump-probe feedline and mechanical oscillator to a thermal bath using the 'in-out' formalism 
of Gardiner and CoUett.'^ Assuming weak system-bath couplings justify making the rotating 
wave approximation (RWA), and furthermore making a Markov approximation for the bath 
dynamics, the following Langevin equations can be derived for the system mode operators 
in the Heisenberg picture: 



ddiiyYl ^ / ^ / \ -J- ^ _|_ 

at n V zmu m 

-Ibmam - i^/2^e'^'^^d^{t) (31) 



10 



and 

dt 



-IpTaT - Zv^e'*-^<(^), (32) 

where 7?,^ is the mechanical oscillator amplitude damping rate due to coupling to the bath, 
7pr is the transmission line mode damping rate due to coupling to the pump-probe line, and 
we have also assumed that the small Duffing coupling and transmission line-mechanical 
oscillator coupling justify applying the RWA to the transmission line mode operator 
terms. The 'in' bath and probe line operators are defined as 

d^{t) = [ c/u;e-'"(*-*o)a,(cj,to), (33) 

v2iT J 



where t > Iq, with the states of the pump-probe line and oscillator bath assigned at to, 
interpreted as the initial time in the past before the measurement commences. For com- 
pleteness, we have also included a classical, external time-dependent force -Fext(t) acting on 
the mechanical oscillator, although we shall not address the force detection sensitivity in the 
present work. 

It will be convenient to work with the Fourier transformed Langevin equations. With 
0{uj) = J^^dcue'^^^Oit), Eqs. (|3T[ and (|32| become 



aM = —. \^/2^e'^''-af{uj)--=l=F,^,{uj) 

m 



00 



+ -^--=^ I du' [aT{uj')a^{uj' — to) + a^{uj')aT{uj' + to)] ^ (34) 



and 



00 — ujt + ilpT I 27r 

fOO /"OO 

X / did' / duj"a^{uj")aT{(^')cLT{uj + to" — uj') 



00 •/ —00 

00 



_^^r^Tm / duj'arp(^u!') [am{uj — uj') + a^{uj' — \ ■ (35) 
v27r 



III. DETECTOR RESPONSE 



The probe line observables are expressed in terms of the 'out' mode operator: 

al'^'it) = /" rfcue-'"(*-*^)ap(c^,ti), (36) 
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where ti > t. The 'out' and 'in' probe operators are related via the following useful identityP^ 

'{t) = -Zv^e-*^--aT(t) + a;(t), (37) 



which allows us to obtain the expectation value of a given observable once aT{t) is deter- 
mined. As illustrative expectation value, we shall consider the variance in the probe line 
reflected current in a given bandwidth 6u! centered about the signal frequency of interest 



Zp J^,-5w/2 27r V (^1 - ^2)Tm 

^\K"M<^M + ar^(^2)arK)), (38) 

where, in addition to the ensemble average, there is also a time average denoted by the 
overbar, with the averaging time taken to be the duration of the measurement Tm, assumed 
much longer than all other timescales associated with the detector dynamics. In particular, 
time averaging is required when -Fext(^) has a deterministic time dependence.^ Expectation 
values of other observables, such as the reflected voltage variance and reflected power are 
simply obtained from Eq. (psj) with appropriate inclusions of the probe line impedance 



From the form of the Kxru coupling term in Eq. (32), we can see that the motion of the 
mechanical resonator modulates the transmission line frequency, and thus a complimentary 
way to transduce displacements besides measuring the current amplitude, is to measure 
the frequency-dependent, relative phase shift between the 'in' pump drive current and 'out' 
probe current using the homodyne detection procedurcl^ While we shall focus on amplitude 
detection, the homodyne method can be straighforwardly addressed and is expected to give 
similar results for the quantum limited detection sensitivity. 



Substituting Eq. (34) into (35), we obtained the following single equation for the trans- 



mission line mode operator a^: 

/oo poo 
did'axii^ — iu')A{iu, Lu') + / duj' B{ijj, uj')aT{u} — uj') 
-oo J — oo 

/oo 
duj" \aT{oj")a^{uj" — uj') + a^{ijj")aT{oj" + tu')] 
-oo 

/•oo nco 

+D{uj) I duj" duj'a^{uj")aT{uj')aT{uj + iu" -uj') + C{uj), (39) 
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where 



B{uj,uj') = 



Tm 



2tX OJ -UJt + IJpT 

1 



Sm{uj') 



UJ' - UJ.rn + hbn 



1 



An U — UJt + ijpT L^' ~ ^rn + 'i'lbii 

St{uj) 



+ 



+ 



-UJ' -UJm- i%n 
1 



-UJ' — UJr, 



C{uj) 



and 



D{uj) 

with mechanical signal operator 



u — UJt -\- i'jpT 

uJrKd 1 
27r UJ — UJt -\- i'jpT ' 



Sm{uj) = ^/2^e"^''-a^{uj) - 



and noise operator 



St{uj) = y^e'f'^-afiu;). 



=Fext(a;) 



(40) 

(41) 
(42) 

(43) 

(44) 
(45) 



For small signal strength, it is assumed that Eq. (39) can be solved as a series expansion 



up to first order in A{uj,uj'), giving the usual linear-response approximation. I.e., 0^(1-1^) ~ 
a'^\uj) + a^\uj), where the noise component a'^\uj) is the solution to Eq. (39) with the 



mechanical signal source term A{uj, uj') set to zero, while the signal component a'^\uj) is the 



part of the solution to Eq. (39) that depends linearly on A{uj,uj'). Thus, from Eq. (37) we 
can express the 'out' probe mode operator as follows: 



a^^'HuJ) 



(46) 



where the first square bracketed term gives the signal contribution to the detector response 
and the second square bracketed term gives the noise contribution. 

As 'in' states, we consider the mechanical oscillator bath to be in a thermal state at tem- 
perature T and the pump line to be in a coherent state centered about the pump frequency 

\{a{u:)})p = exp ^ duja{uj) [a-+(^) - a;(^)]) |0)„ (47) 
where |0)p is the vacuum state and 



a{uj) = -I,\l— 



UJ 



(48) 
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The coherent state coordinate a{uj) is parametrized such that the expectation value of the 
right-propagating 'in' current P^{x,t) with respect to this coherent state has amphtude Jq, 
where 



with Vp = 1/ a/ LpCp the wave propagation velocity in the pump probe line. 

With the pump probe line in a coherent state, we assume that for large drive currents 



Eq. (39) can be approximately solved using a semiclassical, 'mean field' approximation, 
where the quantum fluctuation 6a!^^ (u) in a^^ (uj) {uj) is kept to first order 

only. However, the nonlinear Duffing and transmission line-mechanical oscillator interaction 
terms can give rise to a bistability in the transmission line oscillator dynamics and one must 
be careful when interpreting the results from the mean field approximation when operating 
close to a bifurcation point; large fluctuations can occur in the oscillator amplitude as it 
jumps between the two metastable amplitudes, which are not accounted for in the mean field 
approximation. (See Refs. HQ] and Hi] for respective analyses of the classical and quantum 
oscillator fluctuation dynamics near a bifurcation point). This issue will be further discussed 
in the following sections. 

The solutions to the signal a^r^\uj) and noise terms parallel closely our previous 
calculations, which omitted the Duffing nonlinearityl^ the Duffing {D) term in Eq. (39) 



has a very similar form to the transmission line oscillator coupling (B) term, both involving 
a^a^ operator combinations. We therefore relegate the solution details to the appendix, 
presenting only the essential results in this section. 

The solution to {a^j^\uj)) is sharply peaked about the pump frequency ujp for large Tm and 
so can be approximately expressed as a delta function: {a^\uj)) = x^i^ — ^p)- Substituting 



this expression into Eq. (A3), we obtain for the amplitude x- 

X = c+ [2B{ujp, 0) + D{ujp)] X \x\' , (50) 

with 



7pT ~ iAu y hujp 

where Au = ujp — ujt is the detuning of the pump frequency Up from the transmission line 



resonance frequency ut- Using the expressions for i?(u;p,0) and D{ujp), Eq. (50) can be 
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written as 



Ztc 



27lI^Zp^pT/{hWp), 



where the effective Duffing couphng is defined as 



}C = K,- 



2 -"-Tm- 



(52) 



(53) 



Notice that the interaction between the transmission hne and mechanical oscillator induces 
an additional Duffing nonlinearity (the second term involving Kfrn in ^) in the transmission 



line mode amplitude effective equations of motion (52). However, in contrast with the 



Kd nonlinearity, which can be tuned to have either sign, the former mechanically-induced 
nonlinearity is always negative and thus has a "spring-softening" affect on the transmission 
line mode. Interestingly, by choosing an appropriate compensating "spring hardening" K^, > 
0, the effective Duffing constant /C can in principle be completely suppressed so that the 
next non-vanishing higher order nonlinearity would govern the mode amplitude dynamics. 
Once we have the solution for {a^\uj)), the solutions for the quantum signal a^\uj) and 



quantum noise Sa'ij^\uj) are obtained from Eqs. (A2) and (A4), respectively. These solutions 
can be expressed as follows: 



i^\uj) = ai{uj)A{uj, uo — Up) + a2{uj)A{uj — 2Auj, u — uop) 



(54) 



and 



5aP{u) = (3i{uj)6C{uj) + f32{uj)6C^{2ujp - u), 



(55) 



where the ai{uj) and (3i{uj) functions are defined in Eqs. (|A8[), (A9), (A12), and (A13). 



Substituting Eqs. (54) and (55) into the expression (46) for a°^^{uj) and then evaluating 



the signal component of the detector response we obtairP 
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where the thermal average occupation number for bath mode ui. 

In the hmit of small drive current amplitude Jq 0, we have ai{uj)/c — > 1, a2{uj)/c —>■ 0, 
and we see that the signal spectrum comprises two Lorentzian peaks centered at Up ± uJm- 
The Up + UJm peak corresponds to phase preserving detection, in the sense that 0°*^* gives the 
amplified a'-^ signal, while the cOp—Um peak corresponds to phase conjugating detection, with 
amplifying the a™^ signal.S^l Increasing the drive current amplitude causes the peaks to 



,out 



shift, and the peak widths relative to their height to change, signifying renormalization of 
the mechanical oscillator frequency and damping rate. The noise component of the detector 
response is 



1 



— fku- — 



Up + Au;)2 + 



_^ hus 5u; 



{u-Up + I^uf + 
{u -Up- Au;)2 + 72^ 



|/?2Hl'-Re[AM] + 



{u — Up + Au 




where the integral term involving the (3i{u) functions includes the back reaction noise on 
the mechanical oscillator and the term involving Zp describes the probe line zero-point 
fluctuations added at the output. 

In Sec. |V]we will numerically evaluate Eqs. (56) and (57) and in particular compare the 
detector noise with the minimum noise bound discussed by Cave^^^'^ that follows from the 
Heisenberg uncertainty principle for the detector: 



{[5I°-\us,5u)Y) 

ru^s+5u/2 d^ f ^ 



mm— noise 



hus 5u 



2 2n 



) %T 



X 



X 



7pT 



27r 
27 



Up \u 



Ur. 



bm 



276m, 



u — Up + Au + i'jpT 
u — Up — Au + i^^pT 



{u-Up- UmY + 76m (u^p - - u;m)2 + 72^ 



(58) 
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IV. BISTABILITY CONDITIONS 



We have seen [Hamiltonian (26)] that the current-dependent SQUID effective inductance 
gives rise to a transmission line Duffing type nonlinearity with strength Kd- Furthermore, 
there is a nonlinear couphng with strength KTm between the transmission line and me- 
chanical oscillator. These two nonlinearities correspond respectively to the cubic terms 
proportional to Kd and K"^^ in the mean transmission line coordinate amplitude x equation 



(50). For sufficiently large drive current amplitude Jq and/or coupling strengths KTm-, K^, 



2 

the cubic term x\x\ term in Eq. (52) becomes appreciable, resulting in three real solutions 



over a certain pump frequency range Up. This parameter regime defines the bistable region 
of the detector phase space (the intermediate amplitude solution is unstable and cannot be 
realized in practice). In the following, we determine the conditions on the parameters for 
the bistable region employing the analysis of Ref. [271 

We first express the transmission line mode coordinate in terms of its phase and ampli- 
tude: 

{aPit)) = Me~'^^^'+'t>^'\ 

X = v/2^Me-*% (59) 

where the amplitude M is a positive real constant and recall x is defined through the relation 
{aP{uj)) = x5(a; — ujp). Equation (52) then becomes 



{UT -ujp- t^pT) M + JCutM^ = y2^(&?T)e*^'^^"'"'^'^^ (60) 

where 



Multiplying both sides of Eq. (60) by their complex conjugates and substituting E = M , 



we obtain the following third-order polynomial in E: 



CUt/C LOj.JC^ CU^/C^ 



The bifurcation line in current drive and detuning parameter space that delineates between 
the single solution and bistable solution regions occurs where the susceptibility dE/dujp 
diverges. If we further impose the condition that the transition between the two regions is 
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continuous, i.e. d'^ujp/d'^E = 0, we obtain the bistability onset critical point. From Eq. (62) 
these two requirements can be written 



S/C^E^ + 4{ut - Up)lCE + (cut - ojp? + iIt = 0, 



(63) 
(64) 



Solving these equations simultaneously for E and Auj = ujp — ujt yields the following bista- 
bility onset critical values: 

27pr 



Ebi 



\/3ujt\1C\ 



(65) 



Substituting these critical values into Eq. (|62|) gives 

ypT 



pT/bi 



3V3ujt\IC\ 



(66) 



Finally, using Eq. (61) we obtain the driving current critical amplitude: 



Lbi 



(67) 



'i\/?>UJT\K,\Zp 

Note, the requirement that we operate below the Josephson critical current, /q < /c, gives 
a lower limit on the value of |/C| for which our system can approach the bistability onset. 
The boundary of the bistable region that is given by the vanishing susceptibility equation 



(63) can be expressed in units of the bistability onset critical current 1^ and detuning value 



Acjfoi using Eqs. (67) and (65) to obtainHSl 

3/2 



/ 

hi 



Auj 
2 \Auhi 



1 



Aujbi 
Auj 



± 



AuJbi 
Auj 



1 3/2 ' 



1/2 



(68) 



where the ± roots give the upper and lower boundaries of the bistable region, respectively 
(see Fig. [2]). As mentioned in the preceding section, care must be taken when applying 
our semiclassical, mean field approximations to the detector signal and noise response when 
approaching closely the bifurcation boundary lines. Fluctuation-induced jumps between the 
small and large amplitude solutions of the transmission line mode can occur that are not 
accounted for in the mean field approximation. Nevertheless, in the next two sections we 
shall in some instances evaluate the detector response close the boundaries of the bistability 
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region. For example, we shall see that significant improvements in cooling can be achieved 
provided a way is found to keep the transmission line mode on the low amplitude solution 
branch when operating in the bistable region. 

4F 




da d-i 



I I I 



-3 -2-10 1 
Detuning (lAwbil) 



FIG. 2: Bistable region (shaded) of the cavity-oscillator system for negative, spring softening 
Duffing nonlinearity. The drive current and detuning are expressed in units of the bistability onset 
critical values Ibi and |Aa;bj|. The labelled straight line traces correspond to detection (d) and 
cooling (c) current drive-detuning parameter examples considered in Sees |V] and VI The arrows 
give the direction in which the drive current is varied in order to enter the bistable region on the 
small amplitude branch. 



V. DISPLACEMENT DETECTION 

Assuming that '-/bm ^ IpT, i-e., the unrenormalized mechanical oscillator amplitude damp- 
ing rate is much smaller than the transmission line oscillator amplitude damping rate, then 
the detector spectral noise and response in the mechanical signal bandwidth is approximately 
white over a large range of drive current and detuning parameter space. The mechanical 



signal and noise response spectra are therefore approximately Lorentzian and Eqs. (56) and 



(57) can be parametrized as 



Zp 

signal 
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and 



2mR^ujm ju,,-5u/2 27r (a; - u;p =F Ruj'^mY + (i?776m)^ 



+ ( [5r"*(^, = ^p ± i?^^^, 5a;)] Zp, (70) 

added noise 

where G± is the phase preserving (conjugating) gain (in W-m~^), n^R^Um) is the mechanical 
oscillator's external bath occupation number at its renormalized frequency R^^ujmi R-ylbm 
is the renormalized (i.e., net) mechanical oscillator damping rate, and the detector back 
reaction noise on the oscillator is effectively that of a thermal bath with damping rate 
7back = IbmiRj — 1) and thermal average occupation number n^.^^^- Note, from here on 
we do not consider an external classical force driving the mechanical oscillator; the focus is 



on displacement detection rather than force detection. The added noise term in Eq. (70) 
comprises output noise that is not due to the action of the detector on the mechanical 
oscillator; the added noise is present even when there is no coupling to the mechanical 
oscillator, i.e., when i^Tm = 0. In the absence of the transmission line Duffing nonlinearity, 
the added noise simply consists of the probe line zero-point fluctuations husSu/ (iiTZp). 
However, with the Duffing nonlinearity present, the added noise will be in excess of the 
probe line zero-point fluctuations. 



The convenient Lorentzian parametrization approximations of the mechanical signal (69) 



and noise response spectra (70) that provide the above-described effective thermal descrip- 
tion of the back reaction noise will break down as one approaches arbitrarily closely the 
jump points at the ends of the small or large amplitude transmission line oscillator solution 
branches occuring at the boundaries of the bistable region indicated in Fig. [2j This is a 
consequence of the diverging damping (i.e., ring-down) time of transmission line mode.'^ 



Thus, when numerically solving (56) and (57) to extract the effective thermal properties of 



the detector back reaction, it is important to always check the accuracy of the Lorentzian 
spectrum approximation. 

For sufficiently large gain (i.e., large current drive amplitude), we can neglect the added 
noise contribution and we have for the noise-to-signal response ratio when the mechanical 
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oscillator external bath is at absolute zero [i.e., n{RujU!m) = 0]: 



([^/°"*] ) noise _ .^ 4: . x Tback 

{[OJ^ J /signal 



On the other hand, in the large gain limit the Caves noise lower bound ( 58 ) gives a noise-to- 



signal ratio of one. For large gain, we typically have \2n^^^^^ + 1\ ':$> 1 and thus to approach 
the Caves bound necessarily requires |7back| 'jbm^ 

As an example, we numerically solve for the signal and noise contributions of the detector 



response, Eqs. (56) and (57) respectively, as well as the Caves lower bound on the quantum 
noise (58). We consider Duffing nonlinearities = —3.4 x 10^^ and = (i.e., no 
nonlinearity) . The integrated signal and noise bandwidth is taken to be 6uj = 2i?^7f,m- The 
corresponding parameter values are: probe line impedance Zp = 50 Ohms, transmission line 
mode angular frequency ujt/{2'k) = 5 x 10^ s~^, transmission line mode quality factor Qt = 
ljt/{2jpt) = 300, mechanical frequency ci;m/(27r) = 4 x 10^ s~^, mechanical quality factor 
Qm = uJm/i'^lbm) = 10^, oscillator mass m = 10"^^ kg, Josephson junction critical current 
Jc = 4.5 X 10~^ A, junction capacitance Cj = 10~^^ F, external flux bias $cxt = 0.442 $o, 
and external field in the vicinity of the mechanical resonator B^xt = 0.05 T. These values 
give a zero-point uncertainty A^p = 1.45 x 10~^^ m and transmission line-oscillator coupling 



zp 

Kr-m = 1.1 X 10-5. 

The advantage of using a spring softening nonlinearity, Kd < 0, is clearly evident in 
Fig. [3} where we plot the noise versus response signal under increasing current drive for a 
transmission line both with and without Duffing term driven on resonance, Au = ujp — ujT = 
0. We also plot the response of the nonlinear transmission line for several positively detuned 
values, Auj = ujp — ujT > Termination of the curves indicates the signal value at which the 
damping renormalization i?^ = 0, beyond which the derived solutions become unphysical 
due to the net mechanical damping rate becoming negative and hence the motion unstable 
about the original fixed point. Note that the same criterion, namely > 0, is employed 
throughout the paper in order to ensure stability of the system. Again, the semiclassical, 
mean field approximation is expected to break down in the vicinity of termination points, 
where large fiuctuations in the mechanical oscillator amplitude occur. In Fig. [3} we see that 
with positive detuning, we can further approach the Caves bound. However, this is at the 
expense of reduced gain; depending on one's point of view, large renormalizations of the 
mechanical oscillator damping rate (and frequency) due to detector back action may or may 
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FIG. 3: Detector noise versus signal response at Auj = for harmonic {K^ = 0) transmission line, 
Duffing (Kfi < 0) transmission line (di), and Caves' bound (black-dashed). Noise for the nonlinear 
transmission line is also evaluated for blue detunings: Auj = -|-0.2| Aw;,j| (^2), and -|-0.4| Au;;,j| (ds). 
The labeled curves correspond to the traces in Fig. |2j The dashed, colored lines give Caves' bound 
for the corresponding detuning values. 

not be allowed in detector displacement sensitivity figures of merit, affecting the maximum 
achievable gain as one approaches more closely the Caves bound. 

The trends displayed in Fig. [3] can be partly explained by invoking Fig. [4], which indicates 
qualitatively the force on the mechanical oscillator due to the microwave transmission line 
'ponderomotive radiation pressure' force, both with vanishing and with nonzero Duffing 
nonlinearity and also for 'red' and 'blue' pump frequency detunings. The work done on the 
mechanical oscillator by the radiation pressure force during one period of motion, due to 
the delayed transmission line resonator response, is given by the area enclosed within the 
hysteresis loopP^ESl g^^^j (.g^^ related to the steady-state back action damping rate through 

W 1 

7back = — f^-, (72) 
£/ r 

where W is the work done on the mechanical oscillator, E is the average oscillator energy 
and T is the period of motion. When frequency pulling is taken into account, the usual 
notions of red-detuned {Au < 0) or blue-detuned (Au; > 0) hold only in the weak drive 
limit. We will assume red (blue)-detuned to correspond to drive and detuning values Alo 
where the net work done on the oscillator is negative (positive) as seen in Fig. |4| For a 
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FIG. 4: Cartoon indicating the 'radiation pressure' force exerted on the mechanical oscillator by 
the transmission line mode during one cycle of mechanical motion: (a) the harmonic transmis- 
sion line mode approximation; (b) approaching the onset of bistability. The work done on the 
oscillator is proportional to the area swept out during each cycle, considerably exaggerated here 
for clarity. Positive mechanical damping (red detuning) results on the positive slope side of the 
curves. Negative mechanical damping (blue detuning) results on the negative slope side. A spring 
softening nonlinearity can result in improved cooling for red detuning and improved signal-to-noise 
amplification for blue detuning. 



harmonic transmission line and for low drive powers, the frequency pulling effects can be 



ignored, since the effective Duffing coupling Eq. (53) is proportional to the square of the 
transmission line-mechanical oscillator coupling KTm, which contributes only weakly for 
the considered parameter values. Conversely, the Duffing term causes frequency pulling 
even at low input power and can significantly alter the slope of the response curve. From 



Eq. (72), the decreased slope on the blue detuned side leads to a decrease in the damping 



rate magnitude which, through Eq. (71), leads as demonstrated above to a closer approach 
to the Caves' limit. As mentioned above, benefits in lower noise-to-signal resulting from 
further detuning deep into the blue region are offset by diminished achievable signal gain 
levels. 



Tuning the sign of the Duffing coupling /C (53) to be positive, so that we have a hardening 
spring, results in an increased back action damping rate for blue detuning, and hence a 
corresponding decrease in signal to noise relative to the harmonic transmission line resonator 
detector case. 
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VI. COOLING 



Referring to the parametrizations (69) and (70) of the signal and noise components of 
the detector response, we define the mechanical oscillator's net occupation number through 
the following equation: 

7net(2n^ct + 1) = Ihm [2n{R^UJ^) + 1] + 7back(2nback + 1)' (^3) 

where the net damping rate is 7net = Ibm + 7back = R-yJbm- The oscillator's net occupation 
number is then 

2n„et + l = R-' [2n{R^u^) + 1] + (1 - R-')i2n^^,^ + 1). (74) 

In order to cool a mechanical oscillator to its ground state using detector back action, we 
therefore require a large detector back action damping rate, equivalently large damping rate 
renormalization R^ ^ 1, together with a small detector back action effective occupation 
number n^^^^ <^ 1 . 

Referring to Fig. |4], operating closer to the bistability increases the negative work done 
per cycle on the oscillator by the cavity and hence increases the back action damping rate for 
given current drive. In Fig.|5} we plot the mechanical oscillator damping rate renormalization 
factor R-y, using the same parameter values as in Sec. [V] (e.g.. Duffing coupling Kd = 
—3.4 X 10~^), but with a larger yet still feasible mechanical quality factor = 10^ (which 
we shall adopt throughout this section). We clearly see the enhanced damping as one 



approaches the onset of bistability given by Ibi (67) and Acu^j < (65). 

For the example parameter choices of Sec. |V| we have tOm/lpT ~ 0.5 and thus we are 
operating in the so-called bad cavity limit, where cooling close to the ground state (i.e., 
rinct ^ 1) is not possible )^^ * ^^ * ^^ ! While it is not difficult to achieve the good cavity limit 
> IpT simply by realizing sufficiently large quality factor superconducting microwave 
resonators, together with high frequency mechanical resonators,'^ it is nevertheless worth- 
while to address how nonlinearities can improve on the cooling limits in the bad-cavity case. 
With the fundamental motivation to demonstrate macroscopic quantum behavior, the an- 
ticipated trend is to work with increasingly massive and hence lower frequency oscillators, 
making it progressively more difficult to achieve the good cavity limit. 

In Fig. [6} we plot the dependence of detector's noise effective back action occupation num- 



ber nback on microwave drive current amplitude at the detuning bias Au = —^u^ + 7^^, 
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FIG. 5: Mechanical oscillator damping renormalization factor for detunings both above and 
below the bistable detuning Atjf,j. The amplification region corresponds to negative back action 
damping, i.e., < 1. 



where \Auj\ < \AuJbi\. This is the optimum detuning in the harmonic, transmission Hne oscil- 
lator approximation, i.e., when nonlinear effects are ignored. The noise effective occupation 
number is indicated for both a nonzero {K^ = —3.4 x 10~^) as well as zero {K^ = 0) Duff- 
ing nonlinearity transmission line. We also show for comparison the effective back action 
occupation number when the frequency pulling effects of both the ponderomotive coupling 
KTm and Duffing coupling Kd are neglected. The latter case is obtained by dropping the 



nonlinear microwave mode amplitude term in the mean field equation (50). The sharp rise 



in occupation number and associated sharp drop in damping renormalization at larger cur- 
rent drives is a consequence of crossing over into the amplification region due to negative 
frequency pulling of the cavity response relative to the fixed detuning. The decrease in 
occupation number as / — > is accompanied by weak back action damping, which prevents 
cooling the mechanical oscillator to such occupation numbers. Note that at smaller current 
drives the damping renormalization in the presence of a Duffing nonlinearity peaks above 
the corresponding damping renormalization without the Duffing nonlinearity. This damping 
enhancement can be qualitatively explained with the aid of Fig. |4} In the presence of the 
nonlinearity then, improved cooling can be achieved for smaller current drives. 

According to the above discussion, any improvements in mechanical oscillator cooling 
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FIG. 6: (a) Detector noise effective back action occupation number versus current drive when red- 
detuned at Aw = —^ijj^ + 7pj,, IAlijI < |Aa;f,j|, witli a Duffing nonlinearity (solid line), without 
a Duffing nonlinearity (dashed line), and without both Duffing and ponderomotive nonlinearities 
(dotted line), (b) Oscillator coupling renormalization factor for the corresponding back-action 
occupation number curves. These plots are obtained for the straight line trace labeled c\ in Fig.[2| 



are due solely to enhancements in the detector's back action damping rate for given drive; 
as can be seen from Fig. [6} the absolute minimum attainable detector effective occupation 
number is the same both in the presence and absence of the transmission line resonator 
Duffing nonlinearity. While the effects of enhanced back action damping may be beneficial 
in situations where one is facing constraints on the maximum achievable drive powerj^^ it 
would nevertheless be more significant if reductions in detector effective occupation number 
could similarly be achieved through nonlinear effects. To see how this might be possible, we 
consider detunings corresponding to the pump frequency being to the left and away from the 
cavity resonance, i.e., |Aci;| > |Au;5j|, Au; < 0. For such detunings, the mechanical oscillator 
'sees' a transmission line resonator effective quality factor that is determined by the steeper 
slope on the left side of the response curve (see Fig. |4]). As we drive the transmission line 
resonator towards the lower bistable boundary (see Fig. [2]), the slope of the response curve 
increases sharply and mimics a resonator with larger quality factor, effectively getting closer 
to the good cavity limit and hence resulting in a lower detector occupation number.'^^ESEZl 
Continuing to drive the transmission line resonator into the bistable region, and assuming 
that the resonator can be maintained on the low amplitude, red-detuned solution branch,HSl 
the detector effective occupation number further decreases while the back action damping 
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rate on the mechanical resonator increases (as explained by Fig. |4]). Eventually, the trans- 
mission line resonator becomes unstable at the upper bistable boundary indicated in Fig. |2j 
and the oscillator jumps to the larger amplitude, blue-detuned solution (see Fig. [?]). 
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FIG. 7: Transmission line resonator response curve for Qt = 300 restricted to the small amplitude 
solution branch. The example drive currents are I /lu = 0.8 (green), 0.95 (yellow), 1.15 (red), and 
1.3 (blue). The jump between small (red-detuned) and large (blue-detuned) amplitude solutions 
is indicated by the dotted lines. 



In Fig. |8l we plot the dependence of the detector effective occupation number nback on 



current drive for an example detuning value of Acj = —l^uSl^ + 7pT ~ 1.3Ac(jf,j. Driving 
the nonlinear transmission line resonator towards the upper boundary of the bistable region 
(see Fig. |4]) produces a sharp decrease in detector occupation number, and an occupation 
number value of (2nback+l) ~ 0.55 can be obtained, well below that achievable when ignoring 
frequency pulling effects. The harmonic cavity shows no such decrease in occupation number, 
indicating the qualitatively different quantum dynamical dependencies on and Kj^^a and 
the necessity of the former. We can quantify the effect of frequency pulling by comparing 
with a harmonic transmission line resonator with a quality factor Qt? value chosen so as 
to give the same detector effective occupation number. For the occupation number value 
(272back + 1) ~ 0.55, we have ~ 600, corresponding to cj^/Tpf = 0.95, and therefore 
the mechanical oscillator behaves as if it is coupled to a cavity with double the quality 
factor. This translates into lower net mechanical temperatures as shown in Fig. [9| where we 



give the net oscillator occupation number nnet (74) for various external bath temperatures. 
The combination of nonlinearly-enhanced coupling R-y'-^bm and enhanced transmission line 



27 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1 
Current (A) 

(a) 



2000 



1500 



CE 1000 



500 



X 10 



' I 
' I 
I I 
I I 
I I 
1 1 



1 1.5 
Current (A) 

(b) 



X 10 ' 



FIG. 8: (a) Detector noise effective occupation number versus current drive wlien red detuned at 
Auj = 1.3 AuJhi, corresponding to straiglit line trace c'l in Fig.[2j The Duffing nonlinear transmission 
line resonator occupation number (solid line) rapidly decreases as the resonator is driven towards 
the upper bistable boundary, assuming the resonator can be maintained on the small amplitude 
metastable stable solution branch. In contrast, a harmonic transmission line resonator (dashed 
line) or a cavity with neither Duffing nor ponderomotive nonlinearities in its mean field microwave 
mode equations (dotted line) shows no such decrease in the occupation number, (b) Mechanical 
oscillator damping renormalization factor for the same fixed detuning and drive current range. The 
dashed vertical lines indicate the boundaries of the bistable region for the given transmission line 
resonator parameters. 



effective quality factor can be seen to significantly affect cooling of the mechanical motion, 
even for relatively large external temperatures. 



In the numerical solutions to Eqs. (56) and (57), the Lorentzian parametrizations (69) 



and (70) were found to give good approximations even when the upper bistable boundary 
is approached quite closely. This is a consequence of the wide separation in the relaxation 
rates that determine the line widths of the harmonic transmission line resonator and un- 
renormalized mechanical oscillator modes, i.e., '^bm 7pT- The upper bistable boundary 
has to be approached pretty closely in order for the nonlinear transmission line resonator 
ring-down time to exceed the renormalized mechanical oscillator damping time, resulting in 
the breakdown of the effective thermal description of the detector back reaction. In all of the 
plots shown in this section, the Lorentzian approximation is a good one over the resolvable 
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FIG. 9: (a) Net mechanical occupation number at Aw 



7p-jn for a harmonic transmission 



Hne resonator. External bath temperatures: T = 1 (solid line), 10 (dashed line), 50 (dotted line) 
and 100 (dot-dashed line) mK. (b) Dependence of the net mechanical oscillator occupation number 
on current drive for a Duffing transmission line with detuning Auj = 1.3 AuJbi. The bistable region 
boundaries are indicated by the dashed vertical lines. 



scale of the plots. The actual minimum temperature that can be achieved depends on the 
upper drive threshold where the Lorentzian approximation breaks down, as well as on the 
ability to keep the transmission line resonator on the small amplitude solution branch; the 
latter condition becomes progressively more difficult to satisfy as the upper boundary is ap- 
proached, owing to the increasing probability of noise-induced jumps to the large amplitude 
branch. 

A Duffing transmission line resonator nonlinearity can also produce cooling gains in the 



good cavity limit. In Fig. [10| we consider a transmission line resonator with Qt = 1000, 
giving ijJrallpT = 1-6, and compare the nonlinear transmission line resonator with the har- 
monic resonator approximation at optimal harmonic detuning. Again, by detuning to twice 



the optimal harmonic resonator value, Au = — 2yci;^ + 7p^ ^ LIAuju-, "we see that the 
effective back action occupation number decreases, while the back action damping increases 
as the system is driven towards the upper boundary of the bistable region. Driving a Duff- 
ing transmission line resonator at twice the optimal harmonic detuning can yield a detector 
occupation number (2r;,back + 1) ~ 0.06 just below the upper boundary of the bistable re- 
gion, which is equivalent to an effective harmonic resonator quality factor of Q^? ~ 1400 or 
^m/l^ = 2.2. In comparison, the minimum effective detector occupation number ignoring 
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FIG. 10: (a) Detector noise effective occupation number versus current drive when red-detuned at 
Au) = —yjuj^ + Jprp, \Auj\ < \Aujfyi\, for a Duffing nonlinear (solid line), harmonic (dashed line) 
transmission line, and with the effects of frequency pulling due to both ponderomotive and Duffing 
nonlinearities neglected (dotted line). The vertical dashed lines give the bistable region boundaries 
for the Duffing and harmonic transmission line resonators. This plot is obtained for the straight line 
trace labeled C2 in Fig. [2] (b) Oscillator coupling renormalization factor for optimal harmonic 
detuning, (c) Detector occupation number detuned at twice the harmonic optimum, Auj = 2.2ALi;;,j 
and locked to the lower stable amplitude solution. This plot is obtained for the straight line trace 
labeled C2 in Fig. |2] (d) Corresponding back-action damping rate when driven to the upper bistable 
boundary. 



nonlinear effects is 2nback + l = 0.13. In Fig. 11, we plot the net mechanical occupation num- 
ber for the good cavity transmission line resonator both in the presence and absence of the 
Duffing nonlinearity. Again, we see the strong cooling effects provided by frequency pulling 
of the cavity response. As discussed above, the minimum achievable net occupation number 
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FIG. 11: (a) Net mechanical occupation number at Au = —^Jw^ + ^y^rp for a harmonic transmission 
line resonator. External bath temperatures: T = 1 (solid line), 10 (dashed line), 50 (dotted line) 
and 100(dot-dashed line) mK. (b) Dependence of the net mechanical oscillator occupation number 
on current drive for a Duffing transmission line with detuning Auj = I.^Auju- 

will depend on the threshold drive for which the Lorentzian approximation breaks down, 
as well as on the ability to lock the transmission line resonator onto the small amplitude 
solution branch in the bistable region. 

VII. CONCLUSIONS 

We have provided a quantum analysis of a nonlinear microwave amplifier for displace- 
ment detection and cooling of a mechanical oscillator. The amplifier comprises a microwave 
stripline resonator with embedded dc SQUID. The SQUID gives rise to an effective, Duffing- 
type nonlinearity in the fundamental microwave mode equations, as well as a ponderomotive- 
type coupling between the microwave and fundamental mechanical modes. It was found that 
a spring-softening Duffing nonlinearity enables a closer approach to the standard quantum 
limit for position detection as expressed by the Caves bound, as well as cooling closer to the 
mechanical mode ground state. These findings can be qualitatively explained by considering 
the effects of frequency pulling in the response curve of the transmission line resonator 'pon- 
deromotive force' acting on the mechanical oscillator (see Fig. |4]). With blue detuning, the 
decrease in damping allows for a closer approach to the quantum limit with large amplifier 
gain. Conversely, red detuning towards the bistable point of the force response curve in- 
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creases the back action damping, improving the thermal contact to the detector 'cold load'. 
Furthermore, effectively increasing the cavity quality factor due to the nonlinearity mimics 
the so-called good cavity limit in the harmonic case, allowing cooling closer to the ground 
state. 

The present investigation has by no means exhaustively searched the large parameter 
space of the transmission line resonator-embedded SQUID-mechanical resonator system for 
establishing the optimal displacement detection sensitivity and cooling parameters. Rather, 
our intention has been to point out general trends, using specific parameter values as il- 
lustrative examples. It may be that other choices of parameters (e.g., using a mechanical 
resonator with a smaller quality factor) lead to a closer approach to the standard quantum 
limit, or to cooling closer to the ground state. 

The semiclassical, mean field methods employed in the present work do not take into 
account classical or quantum noise-induced jumps between the small and large amplitude 
metastable solutions that become more likely as the bistability region boundaries are ap- 
proached. Unless ways can be found to keep the transmission line resonator locked onto the 
smaller amplitude solution branch,!^ the predicted effects of nonlinearity-induced cooling 
will be less substantial, as it will be necessary to operate deeper in the bistability region 
to avoid jumps. The driven microwave mode amplitude dynamics in the vicinity of the 
bistable region boundaries is still a relatively unexplored area that requires more sophisti- 
cated theoretical techniques in order to elucidate the fluctuations between the small and 
large amplitude metastable solution branchesJ ^° * ^^ * ^^ * ^° * ^^ * ^^ * ^^ l This will be the subject of a 
future investigation. 
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APPENDIX A: 



In this appendix, we give the derivation of the signal aip\uj) and noise a}^'{uj) terms. 



Suppressing the signal dependent term A{uj^uj') in Eq. (39), we obtain the noise equation 

POO 

a'^\uj) = / duj' B{uj,uj')a^\uj — uj') 



■oo 
oo 



du" 



(0)/ //X (0)+/ // ,^ , (0)+/ //X (0)/ // , ,^ 

/oo /»oo 
duj" / duj'aP^{uj")af{uj')aP{uj + uj" - J) + C{uj). (Al) 
-oo J — oo 



Keeping only terms to first order in A{uo, uj') or a^-* in Eq. (39 ), we obtain the signal equation 

/■oo roc 

a^^\u;) = / du;'aP{LJ-u;')A{LJ,Lj')+ / rfcu'5(w, cu')a?(cu - cj') 



-oo 
oo 



duj" 



+ 



-oo 
oo 



[uj )arp [oj — oj ) + arp' [oj ja^ (cj + uj ) 



(0)- 



duj' B{uj,uj')a^\uj — oj' 



duj" 



a^^\uj")a^^>'^{uj" ~uj') 



+ a^f^{uj")aP{uj" + J) + aP{u")a^?^{uj" - uj') + af^{uj")a^^\uj" + uj') 



+ D{uj) 



duj" I duj' 



(0)+/ II. (0)/ r. (1)/ , // /N 

arp [UJ )arp [UJ jttj, [UJ + UJ — uj ) 



+ a. 



(0)+/ //X (1)/ /■> (0)/ , // /X , (1)+/ //N (0)/ /X (0)/ , // /\ /Ao^ 



Decomposing a^'*(ci;) = {ap\uj)) + 5a^^(c(j) and expanding Eq. (Al) to first order in the 
quantum noise fluctuation 5ap\uj), we obtain the following two equations: 



{aPiuj)) 



duj'B{uj, uj'){ap\uj — uj')) 
duj" 



{a^^\uj")){aP^{uj" -uj')) + (af ^+(cu"))(af (cu" + uj')) 

/OO /*oo 
duj" / duj'{af^{uj")){af{uj')){af{uj + uj" - uj')) 
-oo J — oo 

+ {C{uj)) (A3) 



and 



duj' B{uj,uj')5ap {uj — uj' 



+ 



-oo 
oo 



{af{uj")){a^^^+{uj" -uj')) + {al}!>+{uj")){a^^>{uj" + uj')) 



(0)+/ M\\/JO)f II 



duj' B{uj,uj'){ap\uj — uj')) / duj" 6ap\uj"){ap^^ {uj" — uj')) 
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+ SaP^iu" - uj'){aP{uj")) + 6aP^{uj"){aP{uj" + uj')) 
+ 5aPiuj" + uj'){ai;^^^iuj")) 



+ D{uj) I duj" I duj' 5aP^ {u''){aP {uj')){aP {uo + uo'' — uo')) 

J —oo J —oo 

+ 6aP{uj'){aP^{iu")){aPiu + uj" - uj')) + 5af{uj + uj" - uj'){aP^ {u")){af {uj')) 
+ 5C{u}). 



Mi 



(A4) 



Assuming (a^'*(co')) can be expressed approximately as a delta function, i.e., {a}^'{uj)) 



(0), 



x5{uj — ujp), Eq. (A3) reduces to Eq. (50) for x- The semiclassical approximation to Eq. (A2) 
for aip\uj), with ai^\oj) replaced by {a':^\uj)) = x^{^ — ^p)^ then becomes 

{l-2\xf[B{uj,0) + B{uj,uj-ujp) + D{uj)]}a^^\uj) 

~X^ [2B{uj, uj-ujp) + D{uj)] a^^'^^{2ujp - uj) = xA(cj, uj-ujp). (A5) 

In order to invert and obtain a!^\uj), we require a second, linearly independent equation that 
also depends on a^^^(2c<jp — uj). Such an equation is obtained by making the replacement 
UJ 2ujp — UJ in Eq. (A5) and then taking the adjoint: 



{1 + 2 Ixl^ [B{lj - 2AuJ, 0) + B{uj - 2Auj,uJ~ Up) + D{uj - 2Auj)]] a'^^^{2ujp - uj) 

\2B{uj - 2Au,uj- ujp) + D{uj - 2Auj)] a^^\uj) = -x*A{uj - 2Auj,uj - ujp). (A6) 



Now inverting, we obtain 

a^\uj) = ai{uj)A{uj, UJ — ujp) + a2{uj)A{uj — 2Auj, uj — ujp), 

where 



(AT) 



ai{uj) = V{uj)-^ {1 + 2 Ixl^ [B{uj - 2Auj, 0) + B{uj - 2Auj,uj - Up) + D{uj - 2Auj)]} x, 



a2(^) = -Viij)-' [2Biuj, uj-ujp) + Diij)] IxP X 



(A8) 
(A9) 



and the determinant is given by. 



V{uj) = {1 - 2\x\ [B{uj,0) + B{uj,uj - ujp) + D{uj)]} 

X { 1 + 2 Ixl^ [B{uj - 2Auj, 0) + B{uj - 2Auj, uj - ujp) + D{uj - 2Auj)] } 

+ \xf[2B{uj,uj-ujp) + D{u)][2B{u-2Auj,uj-ujp) + D{uj-2Auj)]. (AlO) 



A similar approach is used to obtain 6a^rp\uj) from Eq. (A4), giving 

SaPiu) = (3i{uj)6C{uj) + /32(cu)5C+(2cUp - uj), 



(All) 
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with 

Piiu) = V{u;)-^ {1 + 2 [B{u - 2Au, 0) + B{u; - 2Auj,uj- ujp) + D{u - 2Au)] \xf} 

(A12) 

and 

= V^uj)-^ [2B{uj, uj-ujp) + D{uj)] x'. (A13) 
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